#execfile("ceorPlot.py")
# TODO:
plt.grid()
plt.rc('font', size=25)
plt.rc('text', usetex=True)
plt.xlabel('$ t $')
plt.ylabel('$ x $')
plt.plot(ttime,xt[0])
plt.plot(ttime,xt[1])
plt.legend(["$x_0$", "$x_1$"], prop={'size':20})
plt.ylim((-1,gas.xlen+1))
outputfilename=prefix+"Fig03-TX_end"+ext
plt.savefig(outputfilename)
plt.clf()
# TODO: 
plt.grid()
plt.rc('text', usetex=True)
plt.xlabel('$ t $')
plt.ylabel('$ v_x $')
plt.plot(ttime,vxt[0])
plt.plot(ttime,vxt[1])
plt.legend(["$v_{x_0}$", "$v_{x_1}$"], prop={'size':20})
plt.ylim((-2,2))
outputfilename=prefix+"Fig04-TVx_end"+ext
plt.savefig(outputfilename)
plt.clf()
# TODO:
plt.grid()
plt.rc('text', usetex=True)
plt.xlabel('$ x $')
plt.ylabel('$ v_x $')
plt.plot(xt[0],vxt[0])
plt.plot(xt[1],vxt[1])
plt.legend(["$p_0$", "$p_1$"], prop={'size':20})
outputfilename=prefix+"Fig05-XVx_end"+ext
plt.savefig(outputfilename)
plt.clf()
# TODO:
plt.grid()
plt.rc('text', usetex=True)
plt.xlabel('$ x_1 - x_0 $')
plt.ylabel('$ E_{pot} $')
plt.plot(xt[1]-xt[0],Epott)
outputfilename=prefix+"Fig06-XEpot_end"+ext
plt.savefig(outputfilename)
plt.clf()
# TODO: 
plt.grid()
plt.rc('text', usetex=True)
plt.xlabel('$ x_1 - x_0 $')
plt.ylabel('$ E_{kin} $')
plt.plot(xt[1]-xt[0],Ekint)
outputfilename=prefix+"Fig07-XEkin_end"+ext
plt.savefig(outputfilename)
plt.clf()
# TODO: 
plt.grid()
plt.rc('text', usetex=True)
plt.xlabel('$ x_1 - x_0 $')
plt.ylabel('$ E_{tot} $')
plt.plot(xt[1]-xt[0],Etott,'o-', markevery=5000)
outputfilename=prefix+"Fig08-XEtot_end"+ext
plt.savefig(outputfilename)
plt.clf()
# TODO: 
plt.grid()
plt.rc('text', usetex=True)
plt.xlabel('$ t $')
plt.plot(ttime,Epott)
plt.plot(ttime,Ekint,'o-', markevery=5000)
plt.plot(ttime,Etott)
plt.legend(["$E_{pot}$","$E_{kin}$","$E_{tot}$"], prop={'size':20})
plt.ylim((-1.0,1.0))
outputfilename=prefix+"Fig09-tEtot_end"+ext
plt.savefig(outputfilename)
plt.show()
plt.clf()


xLJ = np.arange(0.9,2.6,0.01)
yLJ = 4*eps*((sigma/xLJ)**12 - (sigma/xLJ)**6)

